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' . • ■ -' DYNAMIC ANALYSIS -OF 'EVENT HiafbRlES ' . V 

jrhough most sociologipte p-rof ess aii^ interest 'in, social pro^esWas— - V 
how aodial baHavlor and social . systems changfe over ti^e—this'' inteteit . . 
li^iidom^Ta^fiec^^ 



eKitmifle relationships among phenomena at only ana point in time.*" Evan 
i^h^n temporal data are usad/ such as in eKperiments arid in parial'^ studies,- 
sociplogisti' seldom; study . dynamics — the time paths of change*; ^ Xhp;;. 
focus li rat^her on change from one equilibrium level to another , as i ' ; 
measured, for example, by pre- and^ post-differences in an experiment or 
by leyals at successive waves of , a panel* :, "^.i . . " 

: - To soma e3<tent pr^essed interasts an^ 



because data suitable for dyjiami^ analyii Bufe pppprtunitias , 

to collect /such data are often ^passed tiecause'vinvastigaw^ uncertain: 
how to utilise detailed inf drmation ' on change over t:taa. ^ We are BtTUck 
particularly by * the failure .^f social sdientlsti to ^gather and analyze % 
dynamicklly da^a on changes in categdrical variables . , WeL hope , to . stimuia^^ 
interest In the cpllactlon and analysis 6^ what -we tarm event histories. 
An event history recdrds ^ates of events that occur: for iome unit of analysis 

/ ^ -I ' ■ ■ ^ ' 

Examples of. events include changes in categorical variables describing 
ir^dividu^ls—such as marital status, employment status, health status, and 
membership in^ voluntary associatibnsn-as well as , thaae applying, to social ^ ^ 
pollectivities™such as changes in pblitical regimes and outbreaks of 
strikes, riots, and wars. - 




main goals. * Wei wish ito show the value of dynamic /^J^ 

\ ^ J . ^ ■ ^; : ■ : :r i. 

y data for fchfesdcio^ogical stuSy of change in 

■ ^ -r^ ■ - ^ 

[m also wish to describe in some detail a 
\ • . . . / / ^ ; ■ 

king fult MBB/pi. event-history: iata; to answer 



t xntere* sop^al scientists 



V, the .k^ndsl of q\^tl^B| 

Thf (ixpoiiticm * the^ following a^trubtura, Sefction I begins by 



noting the^wide ranga^Mf empirical analyses pdriMtted^ by event^history 



data—includifig^ of 




€^ application of e^^^ing techniques for 



' . analysing ciros^-se^ional^nd pane^-d^ta ,\ 'We argue that the use of 

I continuous-^ time ptochastic models allbws a unified treatment of the out^ 



nd the full use of the iafc 



comes thas^ Variou^/Dtft^^^__a^rMchf and the full use of the ii^f ormation 



in avent^hlstory .data' 



[' '^^^ we take the po 



pogitdan-that event-^history analysis 



h a familiar Tnodal, 'the con^inuous^time Markov model 



be iiiAed to an.'exj^^^ci^ majthemat^ical model .of the social process ' being 
studied * - JWe ^focur^ 

introduced intb spci^ogy by^blenfan (1964) ^ We suimnariza properties of 
the model ■ and its iinplicatidns for. various outcoifias , wl also describe 

■ ■ V . ^ f, . / , y\ ' ' 

extensions of this model/co deal WiGh^population hetalogeneity and time- 
dependence.; these eKtensions makS the. model much more useful ^'f or sociological 
researchp . Sec,tion oo^ludes wlth /^n outline of a maxiii^m. likelihood pro— , 
cedure for usirm..-?event^history ..data 'po estimate and test the extended model. 



, _ IL 

Section II givfis an empirical llluacration. This involves an elaboration 
of tfie analysis presented in, Hannan , Tuma , ' and Groeneveld ('1977) of the 
effects* of* an income maintenance experiment on e'hange in marital status* 

^ . . . . ^ . ... ..^ . . . . 

■^To illustrat.e the extension dealing with population heterogeneity, we 
report estimates of the efifects of experimental Weatments on rites' of 
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■J 



marifiagej tittarltal dlissplu^ and fexperimantal^ atirition. fo illusttatf ^^ 
the exCenp Ibn dealtng \^ith. time-y^pendence j^^^w^ estiraatas of hiw 

eatmeat effects vary oyer the exper imen'talj* peridd . ^ We conclude ' by 

about' the process than; oth^r procedures ' more familiar to sociologists 
but also fits the. data; well^-at least as well as these alternative 
procedures.. : ■ ^ " . \ 

Fiftallyi in Section III we contrast .dynamic analysis of event -^history 



data with other kinds of analysis. In par^#cular, we compare it to cross- 

^ ' ■ -\ ^ " ■ ■ . ' i ■ ■ . . ^ ■■ 

< sectional analysis (assuming the process is- in equilibrium), evant^count 

^ & ^ ' . ^ ■ ■ ■■ . 

analysis j and panel analysis. We find that when, event-hi$tory data can 
be obtained, the strategy and_ procedures that we advocate have -substantial 
' advantages ove hese thr^=^ alternatives , ' ^" '■\ , / ■ ^. 

Obviously a paper of this sort builds on a'\,nut^ber of intellectual ; 

' ^ ^ ■ '^'^^'^^ \ ' - 

"V^-feMditions and on contributions of previous auth8r%l Our debt to 

classical mathematical sociology is obvious in hxix d^iscussion of this 

. . ■ ^» ^ i 

. -^^ y ' ; • ■ 

model aT\d its- implicacions . Readers familiar with^-the t^^al^mgdeling 
(or structural equations) tradition 'should discern its influence on out 
• approach to modeling causall structure and t^me-dependence . In a cbmpaoion 
paper (Tuma and Hanna^ ^ 1978) , we. ha^^e sketched the genealogy of our 
statistical apprQach and have briefly reviewed some alternative estimation, 
strategies* , ' ' * j 



CHOICE OF A MODEL 



^k3r^i;tm'T^ 

analyse event historiej, we have even more options than u^iuai. This is 
eijly' Irf^thi various dependent variables that can be selected from 



hiitories. *foT example, consider how differenc investigators interested 



J) 



seen 

^ in ^ni4l:ita^stab^lity might u cilice . Marital histories to itudy cfianges in 
marital status*^ "Sorte-^VouM study the probability of at les^st one 
marita^^issolutiori' within ^me fixed period-. Others mi|ht study the 
average length of marriage*. Others might choose to study the number of 
'marital events during a fixed pariOT, And, some might analyze change in 

' ^ ^ ■ ' ■ ( J : ^ . ■ ^ ■ • . . . . , 

£he proportion %f^a population that is married. " 

^ " ■ k — ^ ' ^ ■ . ^ . ■ ' 

, Whi.t/are we' to make of the richness of the informatiion c6ntained^in 
'eve.nt histories? We certainly have no objection to different analysts ^ ^ 
5 focusing on different, aspects of 4 problem, fiowever^v thd various dutcomes |- 

^ " ■ • > ' . - >- ^ ' ' ■ ^<^-- 

listed above are hoc independent. Moreover-s no^'ione of t'hem 'conveys all of 

■ ^ ^ - ■ -1^- • ^ 

the inforfcation contained in an event history . ^ We bel^4.eve that much 

; • \.\ ' ^ ' ^fc . ^ ^ % ■ - ; 

is gained from/designing analyses that-ytilize all of the" ^ 
information in Un event history, /th^t is, which use inf o'^rma^oji on the ^^^4 



nuiefoer, seqiience and timing ^f events. 



. * Our solution to this pjoblem — not necessarily ^fihe pnly one "rests upon \ 
' ^ the recognition that w% can infer each of the ^fependent variables listed^ earlier \ \ 
¥ if we can describe each person's stptus at tevery 'moment in tj^me* More gerterally^ '-^ 
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we can uni^ the analysis of various' outcomes that caft^ba aktracted from 



variable^at every Mment^n .time Th strategy has' the advantage o£ pareimQny--* 
'one model replaced several. It also has the advantage of dealing with the 
intier dependence among .the^ vatfiou's ^outtome&* Of courses ^® expect to pay a 
pripe 'for ^uch analytic poweor- 'In partibulax, ^ wa must 'make some restrictive - 
assumptions about the process generatirfg the^ev^^^ts. l^ecause ^f th^ae 
assumptions, our model may . fit sample data less well""£or any particular 
outcome at tfny arbitrary time than a model designed spec^if ically to^ccount ^ f 
for^ t'hat . buf come , Ajfrter we have desctibed our model and estimation^ procedure, 



we^ire^utn to this is/sve as part of an^eKtended empirlcfl illustration of our 



! methodoidgy* 

A Markov Model .of Events 



i^W fii1.lv_ tj^ m1!.^^ ^f^T^.ulate an* explicit 



' To use event his'torieV fully, we must fo 
model Q'fr^' the process generating Ihe event^^Because' the model^ ^ ^ - j; 

must account for states of ' ^ategoi;ical variables' at ^ery moment in time,^ 

■ - _/ ; ^ : \ : ' . , ■ ; - - : 

and because* the events can occur at - any tinte. We choose con- ^ 
J ^ tinuous time models* Exactly wihat' type of raoiel should be ^sed depends 
. ^ M ) on thg substantive .application ; variety of con tinuous'-- time /mo'dels could 



f 



'be employed with the general rriethbdology w^^ptopose. Her,e we c^centrate 



# 



m one type af model* a f inite^state, continuous- time' Markov model. 

. ' / ... , r i - . 3 

.V^begirf by briei4y gtating the formal assumptd/bns of ^thl^ model. 



1 ' 



( - 



* r • - 
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. - , ..... . ^ ■ ■■ . ■ ■ ■ 



^ Let Y(t) be a tandop.vwiable denocing the state of .the categorical variable 

' ' ' ' - ' \ ^ ' ' * ' * = 

occupied by a member of a ^homo^^eneotts popqlation at time t. The get of ail 

posgl ^ie_val u es of Y(t) is tailed /the ^tate- gyc e o f t Ke p rocess;' this is ^ 

issumed to h¥ finite. Let p.i Cust) represent the probability that somaone in 

Sjtata j at time u is in state k at tima t, ti*- ' . 

, p j^Cu,e|( - arob[t(t) s k|Y(u) = j] ; ^ CD 

' , ' : ■ ' ' ' 

and let P(u,t) denote tha*^matrlx» of-'SEhese transition probabilities." 

^ P(u^s c) ■ ^IPjj^Ou, t)} * For example^ suppose that j k = married; than 

p22Wjt) stands for the probability that those iriarrild ai time u.ara^also 

rtamed^ at= tnjna t, ^ ,i ^* ^ 

"^"^ . ^ . ' ■ : 

' Next we dafina r , (t) as the' Ins tantaneous rate of transition from 

jk ' . ^ « 

f,ata i to stace^k at time t. Ifce transition rate r\, (t) is the limit, 
as 4t ap^roachas zero, of the^^^probability . of a^ change .^from j to k between 
t and t+Atj per unit of ^ima- ' - ^ K j . 

■ " * I ^ ■ . ■ ' ^ ■ - ^ ; 



r (t) = , lim ^ j#k . ' (2) 



The rate of leaving state j at tima t, r.(t)*, is ^s follows 



By * the Markov assumption, *the Chapmari-Kolmogorov identity holds: 



\ > ■ P(v,t) ^ BCv,u)P(ijp^) v^'u^t . (4) 
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Given this,, ana Che usual addltional^assumptlons 



continuJ^ 



; ) 



dP(u.c) 
dt 



PCu,t)R(t) , 



(5) 



where R if .a matrix in vhich the j-fe of f-^d^aional element is the tr^sition 

th ' 

^jk^^^* the j ^diagonal element is the negative of the rate of, 

• - J ^ ' ' _ ■ ' : ^ 

leavings -rj(t) . In the time-independent case, i.e,j R(tO ^ equacioh 

4(5) has Che solution: * ' / ' 



(5) 



where :a ^matrlx raised . to the zero power is defined to equal the identity 
maCrlx I (so ^ I) * _ # . 



The^irst task in using suah a model is specification of the ^tate 
space^ i.e., the exhaustive and mutially eKclusive set of^screte values 
of Y(t)/ j^For example, in our analysis of the 'effects of incoTne maintenance 
treatments on marital stabiiity (see Section II below) j the relevant states 
of Y(t) area (-not married) , 2 married) , 'aad 3 ( -latlrited) , State 3 
includes^ those who refuse to participate in the Stndy /cannot be located, 
die or:, emigrate^ outside the continental U.S. "ittrited" is an absorbing 
state; it cannot bereft, Thur R has the formi ' . : ^ ■ " 



.R ^ 



^r. 



^12 ^13 



^21 "-2 



^23 



(7) ^ ^ 



AO 



Although the elements P(u,t) cahriot be written . as explicit functions of 
the transition rates , for a general matrl^ R, . this can done %en R has, the . 
form In (7)* In particular* it can be shown that P(u,tO has the elements i \, 



, - r X (t-u) 



(8)' 



X^(t-u) . X^Ct-u) 

e ' - a 



X^Ct-u) XjCt-u) 




(10) 



p^j^Cu.t) = p^^'CujC) =,0 



(11) 



P33(u,t) 



= 1 



(12) 



where j=L or. 2, k=3-j , X^^\^^ and 



/2 



(13) 



^2 — , 



•^+r2=>/(r^=r2)2+4r^2^ 



21 



(14) 



The explicit incluaiofr of attrition as a state is an imporcaiiC feature of our 
applicatlQn of. a Markov ^model. All too frequently 'invastigators with temporal data 
only 4naly^ze the fraction of the original sampla'that has not attrited by 

the J ime of some later measuramenc . This procedure implicitly , 

assumes that attritionJias no impatt on the assessment of the effecGs of ' 



o 
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Gausal variables, BKpllcitly inalMing attrition as a .state, we avoid 



this -highly dubious aisumption. HowevOTj^it is still necessary to assume 
^Chatr^Ee^^¥Se^mo3^r^app^^^ 

do npt* This assumption could b& incqr^ect^ but cannot "be 'tested 
without a follo^up study of marital status changes of thosa who /. , 
attrite. \ . ^ . 

The inclusion^ of attrition as a state complicates constdarably both 
the derivation of equations and our discussion. For this reasons parts 
of our discussiou are based on the twors tat e model i'n which there is no 



^ 1 



attritions 



R ^ 



=r. 



12 



Than we obtain 



p (Us t) ' ^ 



Ik 



1 ^ p . . (u, t) ^ 



12 21/ 



1 - e 



t-u)j 



for- jjk ^ Is 2 and 



ications of the Model 



We^ claimed earlier that use of such a eontinuous-time stochaatic model 
lets us' derive implications 'about a- variety of observable variables, .In 
this sectiari show thiis by discussing a number of well-known results 
about Markov mbdels. " ^ ^ \ " ' : * ' . ^ 

Time between Transitions .- It il-w^ery kno^Cthat for a^'^^Jcov model 
the length of time between transitions has an exponential distribution whose 
paramcti^r 'depends on the transition rates (sees e.g. ^ Breiman, 1969.,^ Ch. 7). 



10 



In particular, let F, (t |u) ■ repre,senc the probabiiity of a transition from 

state j before time given state j is occupied at time u. F.(t|u) is 

usual-ly called the cumulative probability distributioTL function, ^or simply,' 

the probability di-stributidn function. We can show that for' modelfe iti 

, which R(t) ^ R,^ . , * ' 

* - -(^t--u)r, ^ ' 

F. (t ju) ^ l^e . ^ ^ * 

Later we use the definition that , ^ 

■ 

(t |u) ^ 1 - (t|u) 

G,(t|u). is often called the Survivor function becausa it gives the ' 
probability that a unit in state j at time u remains, (or survives) in 
state Ujntil time t. For models in which R(t) = R, 



3 



-^(t-u)r . -^(t- 
J 

u) ^ e - e 



This says, for example, that the probability chat a marriage existini 
at time u survives until a later time t declines axponentially as the 
length of the interval (t-u) increases. This monoconic decline occurs 
because it becomes increasingly likely as time t increases that either 
the marriage breaks up or the couple drops from the experiment. 
We also use the probability density function: 



. dF. (tlu) - dG, (c|u) 

f (t|u) S J . j • 

J ■ dc dc 



9 
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Jt can be shown that in general 

. , ' fj(t|u) = r.(t) GjCtHu). . ■ _ ' ,(21) 

For tha.parci'cuiar^ase in which kit) = R^" fhao la, the-^^ffisitito^ rates are 
constant over time , • ' _ 

-Cc-u)r, 

f J (t |u) * r^ a J ■ ■ (22) 

The probability of a change from scate j to some other state between t and 

(t'+ dt) is apprpximately equal to f ,(t|u)dt. Equation (22) sho^s that 

this probabiliCx initially (I.e. at (t-u) - 0) equals r and dacllnes ' 

. ' * J ^ ■ 

exponentially as the length of -the interval (t^-u) increases. In other words, 

the probability of leaving a state varies over time even when transit|.Qn rates 

.are constant. This is one of the main advantages of modeling social proceSjSes . 

in terms of cransition rates rather than probabilities of change, 

Xhe average duration of state pccupancies (e.g., expected duration of 

marriages and of the intervals between marriages) when R(t) ^ R is easily^ 

shown to be ' . ' 

* E(t-u) f 1/r . , (23) 

where u and t denote the times of entry and exit, respectively, from state j,^ . 

Given that a. change occurs at time t, n^x^ Ib "the conditional probabi= 
lity that k is. the destination. Thus, with fhis moder, we can account for 
both the probability of at least one change within any specified time 
interyal. and the conditional probability that a change involvest"TOyemenc . 
into some specified state. ' ^ 



12 



State probabilities .' We can also use the model to find the unconditional 
probability of being in any speciflea state at atiy point in time (e^g.^, the 
probability of being married at any specified time). Let P^Ct) denote the 
unconditional probability of being in state k at time t^ and £(t) represent a 
row vector giving the probability distribution among the states at time t. 
It is easy to see thst 



(24) 



or 



,E(t) ='£(u) P(u,t) 



(25) 



Usually we are interested in £(t) when £(u) = £(0), the distribution at 
the start of the procesS;; .; 



For simplicity we begin by considering the two-Btate moder, e.g., 
for Che study of marital stability when there is no attrition. Inserting 
equation (16) into (24), we obtain: 



'21 



'12^^21 



1 - a 



-Cr +r,^)t 



(26) 



?,(t) ^ Pt(0) 



■"^-12^^21^^ 



12 



^12^^: 



1 - 



T--12*-21^^ 



(27) 
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As time t becomes very large, Cc) approaches a stsady^otate (or cquilibrimn) 

value 1 p^(^) "'^21^ ^^12^^21^ ' ^2^"^^ " ^12^ ^^12^^21^ PjCt)'ls^a weightad averagi 
of the initial proportion in state j and the s€eady=state' probability in state j. 
The weight given to Pj (0) declines exponentially over time, while the weight given 
to the steady-state proportion in state j .increases over time until it reaches onf . 



For the three-state model, e.g,, for the study of marltai'^stability when there 
is attrition, we find (by substituting equations (8) through (12) into (24)) that; 



(r2+X2)p^(0) + r2^(l-Pj_(0)) 



- e 



^1 = 



{(r^+,\pp^(0) + r,^.(l=p^'(0))} 



(t) = 



A,t: 



{(r^+X-)(l-p^(0)) + ^^2^i(0)}^ 



' B {(r^+X^)Cl-p^(0)) + r,^p, (0)} 



12^1 



/(X2-X^) 



(29) 



(t) = 1 + 



X,t 



e {X^ + ^13?!^°^ + r23(l-p^(0))} 



X,t 



/(X2-X^) 



(30) 



where; X^ and X^' ara given hy (13) and (14), respectively. It is important 
to emphasize that r^ = r^^ + r^^, = + r^^ and that X^ and X2 

depend on and r.. .' Thus, each of the above equations depends on Che 



actrit ion ^ates , r^^^ ano ^73' 
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Expacted number of events* Finally, we can* derive relationahips between the 
fundamental parameters of the model .and the expected number of events that occur in 
any Interval of timei This, requires extension of the ordinary Poisson model of the 

ttumber of events in a time interval^ since the Poisson model assumes . that only 

one kind of event can occiff. Let n.^COst) denote the number of changes 

to state k between times 0 and t when state j is ^occupied at time 0, ^ ^ 

0 ^ t* For example^ n^^COjt) stands for ithe- number of marriages formed 

before t tor a person who da married at time 0, It can be shown that 

for^ the two^state model (when there is no attrition) \ . ' ^ 



^12^21^ ^12^21 

■ un^CT*— -,2- 



J.- -J.., ^ 12 21 



e ^ 1 



^12^21= 



(r +r ) ; 

12 21 



e ' - 1 



(32) 



where J,k = 1,2 and j#k. For the three-estate model we find that: 



X t A. t 

(1-e - ) (1-e ) 



X, 



A, 



(33) 



A,r A c . 

Jrj^+Aj) (1-e ) (rj^4.A^)(l-e - ) 



(34) 



j,k = 1,2 and jj^k. 
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. \ : - : 

The point we ^wlsh Wo emphasize is 'that these diverse observable measures' 
on event' histories are all explicit functions. of the rates that, define the 
Markbv model* It is not necessary that'each of" these relationships be 
estimated separately. Any one of tb^m' provides an approach to estimation. 
Once the .rates have *be^en estimated, the^remaining quantities can be calcu/lated 
Using these functiq;ns. The 'Empirical analysis in Section II below uses these 
in #valuating the fit of our model. - , 

Two EKtenjlons of the ypdel 

Markov models have been used quite ofcen by sociologists in the past 
'twenty years or so. In the simple form discussed above, these models have> 
invariably failed to fit the data. Sociologists^ efforts to increase realism 
and improve fit have concentrated on modif icatio.n of two assumptions of the 
simple Markov rriodel^ ^. . , ' . ^ 

constant transition races. Instead, fdUomng Coleman, (1964^ 1973) and Tuma 
(1976), we may assume chat Che same, constanc rates only" govern the be^havior 
of .units wiCh ■iaenciaal- values on a sec of obsar^le, exogenous , variables . 

In other woras , we may es'tablish relationships becwe^n observable variables, 

■/ . ■ ' - . . 

denoted by X," and the rates of dhange in caceiorical variables.. The mam ^ • 

socioloeical interesC in mosc analyses of event histories lies m cesting_ 
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hypotheses concerning the e'ff^ts of aSogenoua variables on rates.' For,* 
example, in analyses' marita^ histories produced by the inaome maintenance 
;^xpetiments, we concentrate on ^testing hypothas.es about tffaGts of experimental 
treatments on rates of ma'rijtal dissolution artd reflarridge. / 

To introduce j^such -canisal relationships, we must state pn 'explicit^ . 

- . ' _ ■ ! ' ' ' ' " . o ' . ' " ■ 

dependence of the anobsfe^vable ^atjss on^ the observable variables. ' 

In our empir^tal analysis j we assume a Log-lineat Y^i^ciQn|hip between each , ^ 

transition rate and'X: - 



^ri r 



X 1^ for all j and kj j^k ^ (36) 



where 9^,. represehts a vector' of ^parameters to be estimated. . The log^ ^ 

" J *^ " . . . ^ . 

_^ . _ _ - 

linear relationshJ.p constrains r^^ to , be .positive for each individual, . 

whatever the value of in accord with equation (2); we. find that it.usually 
fits data better than a linear relationship. We also assume^ that r ^i-S 

finite for air j'and k and for all Individuals. * ' \ 

Tlm6 -dependence , we have assumed to this point chat the rates 
are constant over time. Sociologists are accustopied^ to thinking about 
population heterogeneity (i.-e-,- causal relationships ) ^ but have not devested 
much attention to^ issues involving time=-dependenc^ of effects,, l,ie,,,the 

mariner in which causal relationships change over tipe. However, sometimes theory 

^ . \ '-'^ I , 

\ indicates.' that rates of change are some specific function of time. For example. 



, ^ -J ■ 



Sorensan .(19755 has argued that rates ai_laaving a job are exponentially f 
declininl functiorfs labor'' fbrce experietice ^ • while'^ma' (1976) has ' ' ^ 
specified model in which rates of job leaving are secdnd'^ojjder ^ * 



polynomials in dur^atio.n i^^xhe j^.b* Such paramet?Tc forms of time-dep%ndlnce 
.may' be'' acc^^^acad lr\ ^chj^ strategy discuss without difficulty. 

, ' Sometimes^ ^o not have any ^ a priori hypothesis concerning time- 
dependence. In an ejcperiment , we migKt exp.ect that transition rates during 
an initial adjus'tment' period 'differ from rates later on, but' be unsure whether 
initiff^rates should be larger or smaller than later rates* ■ Tm these circum- 
stances we can define a set of time pe^iodfe and allow t^e ef¥ects of experi- 
mental variables to vary^fteely frM period ^£o period^ while remaining constant 

within each period, In^ thi% situat.i6n t^e -time periods may be arbitrary^ a 

•> ^ ■ ■ ^ ^ . / ' V 

disadvantage of this approach/ Hqwevers ,in other sitiMtions we .may have some 

^idea about points on the time axis when ratefe'may change In some way. Fo^^^ 

example, rates of collectiv^ violence may change when there is a chang# in ^ 

the political regime or during periods of warfare, . , . 

■ ' ' ■ ' 

"■iMaChematically stated, a model in which rataa vary from one cime-period 

■ A V ^ ' ■ ; 

to ariotlher is: 

' 1 * , • , ■ ' ■ ' 

; jkp -jkp-p . p-1 - p . (37) 

■ ■ - , 

where t_ is the last moment in peMod P. Tq = 0. is the- starting time, Z_ is • 
a vector pj exo|enpus va'Tiables affeccttng the rate in period' p, and is 
a vector of paramacers giving effects' on|rjj^p in period p. We illustrate ' 
the use of a model ilmirar to this one belbw. 
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"^.phe .linea^^^rugtural, aqulfc-lo|is modajs^ith .which socioWgi^s have 

.becom^ quite familiai 'reefeijt years have-Tfl%3e eoTranon the 'view chat the ^ t * 

« ^ . astlmation equatlgn i4/iaentical Bu^r^fthis ±i not necessafilT . 

- ^ .the' case,^" The equations rsrating-obsrfrvab^ outcomes (e.g, ^(^averg^W duration, ' 
/^^^ ^ ^ . ^ ^ ; ^ ' ' ^ ^ ■ ^ . 

^number of events, ^etc*) to the rates are^ll .possibla candidates for^ 

estimation equations, but they are not equally promising because of their ^ 

■ ^ ^ compl^ity, ' * J 

' ^ . / Building upon .the work of Mrthoiomfew (1957) arid Albert (1962), we \^ 

form maximuirt likelihood^ estimatorfe of rat^ (and of%ausal effects on ' ' 

rates) using data on the da€bs and kind of e-^ents. Ihis ap^oach offers a 

^ - _ . " ' ' ■ ■■ ' ^ ■ , 

lumber of advantages in the present circumstances. First maximum likelihood " ^ 



^tlmators have good^ la£^ge\ sample properties under fairl^ general conditions* 
Sepond, maximum likelihood estimators retain their good properties under any 



mpnotonic tranBtormation* Thus, one can use maximum likelihood estimarers of . 
rates do form maximum likelihaod estimacors for expected durations and * '* 

other monotonic functions of f^tes . Third, maximum likelihood procedures ^ 
permit a satisfactory solution of what is called the censoring problem. Data in 
which ttte observation period is too short ^to -record a change for every cas« 
are said to be censored , and errors of inference Ire likely if appropriate 
measures for dealing with censoring are\nbt adopted (Tuma and 'Hannan, 1978), 

Assuming independent observations' on the N different cases (units)^eingf = 
analyzed,^ we can write a likelihood function that ,uses^ all information in 
©vent histories-'^namely , the number, timing, and sequence of events. This 
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inf6rmatiDn can be represented by the^ following kinds of varia-tlesf 

occurs in^^he 
and 



dummy variable chat equals unity- ^if .case i's m - event occurs 

\ ' • ° • *• ' 't-'h /• 

observation perJLod), t '.' (the cline of case i's m event if w =\ 
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the 'time that casa^^i Is^ast > obs^Jre 



Chat is unity if casff i's m' ' event is bbservey;and consists of a changa to 
state j-and acherwifee . equals zero). We 'define t.: to be the start of the 



d if ^j^'^'O)^ and v^j^ (a dummy variable 



' ^. o^^ervation|^periQd on case i and y^^^ to equal unity ii state-.j' is occupied 
'^y 1^ at t^^ and otherwise ^zero-. Note that I , w. . equals . the 'toil number 

,of evants observed co.>occur to 'case"'!, while^the t 's are -the times that 

' ^ \. ■ " " 

.these ei(ertts occur ind che Y indicate -€lte sequence of events. 

'iV^.. simplicity we begin with the l^kelilfiod equation that uies infortna- 



t^on on the first event only; 

f.N n 

= IT n 

i=l j=l 



[ 



(1-w,)Vq. 



n 

k=l 



r., W/r.(X) 



(38) 



wher^ N is the number of case=! r, i= 

_ . - r case^, n is the number of states, and the subscripc 

i on variables (namely, w^^ . . , ^ ^^^^ ^^^^^^^^^^ 

clarity. The first t^erm [in square brackets], ,ha survivor function, gives ' 
.he probability .hat thefi.st event has not occurred by ti.e t^ (se. equation 
(IB)). Since any number raised to the .e;^^^ ^^^^^ 
term differs from unity only for those cases that have not had a firste.ent 

(w^ - 0) and are in state 1 at time t fv = n , Th 

() ^ Oj " • The second term is the 

probability density that the first evant occurs ^t tW r Y • .n 

-111= occ^urs^c time t^ (see equation (20)) . 
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This term differs from, u^ity only far thd^^cases that haye had 
M'^kn.t (w. ^ 1) and a.re in state 



20 

at leasc^ onm 



j at time tg '(vq^ ^ 1). khB third ^te^^s 



probabil4.ty that the fifst event consist^ of a move to state k, given that state 



^J^^ii^ lef t ^(see equation \;36)). This terms differs from unity only for those cases 
■ hav^iiad at least one event {^^ ^ 1) ^ that are in state j at time Cq (Vq^ ^ 1)_ 

and whose first event consists of a move to state (v z^l^) * 

* ^ V ^ Ik 

' ■' Equation (38) q#n be written more simply by sXibstituting (2) into (38) 
' and, by Chen collectlrig^arms : , - - ' . \ . 



th 



i=i j=r 



oj 



n 

n 

k=I 



W V V 

-I 1 OJ Ik 



WTien R(C) = R, so that (19) holds., we can simplify this sti]il more; 



N n I n 

n n ] .n 

1=1 j=i k-r 



(39) 



OJ 



n 

a 

k=l 



Voj^lk 



^ n n 

1-..= n n 

J-l k=l 



n ■ 
1=1 I- 



-(t^-tQ)r.j^(X) 



'Oj 



^5l^0j^lki 



Note that the term in curly brackets in (40) above does not - depend on 
where j ' # j and k ' ^ k> Therefore, the maximum likelihood estimates 

of a given r,, -^or of the effects of causal variables on r.,as in (36) " 

jk ...... jk 

/can be obtained ty maximizing the term in curly, brackets for that particular 
j and k. This fact is quite important in practice. It means that we can 
estimate selected transition rates without having data. on all kinds of 

. transitions. To escimace r^^,, we, need data only on the' times of all first 



(40) 
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( avents (t, ) for cases originally in state j, (v = 1) and*on the outcome 

of this event (the value' of ^j^^) * This, means we can concfntrate data 

collection and analysis on urritg originally in states of particualar theoretical 

interest and can ignore events occuring to- other units ^"^ unless we wish to 

model or predict the overall evolution of the entire process* 

t 

Equations (38) through (40) can be generalized in a straightf brward 

r 

way to th© situation in which data on a_l,l oSs^rved events (and not just the 
• first) are us^d in the estimation procedure. We give the equations corresponding 
to (38) and (40) only^ leaving ,che intermediary aquation to be supplied 
by the interested reader: 



N 38 n 

n n n 

i-l m^l j-i 



G (t_/: . t X) 



n (l=w )v ^ . 



(t I t , , X)^ 



w V . ^ 
m 1 5 j 



r,, (X)/r.(X) 



w V =, ,v 
m m-^ljj mk 



(41) 



^'hen R(t;)^^ R, this simplifies . 



n n 



j-1 k-1 i-r 



-(t ^3r., (X) 



\ 



^ 1 * n w V -f . V , I 

5 J ^ I ' T m-1 , ] mkv 



Maximum likelihood estimators of the parametera are found by maximizing ^£ 
■^ot its logaritta).^ , The optimal asymptotic properties of maximum likelihood 
estimators .(consistency , asymptotic normality) are well-known (e,g,, see Dhrymes, 
1910). There is also evidence (Tuma and Hannan, 1978) that the properties 
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' ... ' ' ^ 

of the eBtimators obcainad ' f rom (40) remain quica good even in small samples 
and with a higli degree of censoring (i.e., the mean of number of events is- low) 

= A likelihpod *ra^io Ce#t can be used to test nested jnodels* Let^ be^ - 
the likelihood ^for^a model Q with q+s estimated parameters and ^ be the 
likelihood for the nested model that has .s parameters constrained (usually 
to equal zero) and q ^parameters' estimated . The likelihood ratio X ±g defined'' 
to equal max iXj/ma^ (^^. Asymptotipally ^7 In X has a chi-square 
distribution with q degrees of freedom, permitting us to test the fit of 
the model Lu relative to the model Furtherinore, it is possible to perform 

tests bn 'the' coef f icients of individual variables using the estimated 
covariance matrix of the parameters (the inverse of the matrix of ' second ■ 
derivativefe of the natural logarithm of ^ with respect to the parameters). 
The square root of the variance of a parameter gives ^its estimated standard 
error, which can then be used in a standard fashion to caLlculate a 

t^statistic (or F--ratio'* that tests whether the parameter differs from its 

valuejf" n the null h^/po thesis. 
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v^II ^ M APPLICATION . ^ / ' ' ■ 

ji* ^ \ ^ ' - 

To* illustrate the\_models and methods discussed abov,p we apply 'them^ 
to data'on white women during the first two years of the Seattle and Denvet ' 
Incom"e Maintenance Experiments {SIME/DIME) . We combine data fro'm Seattl^e ^ . ^= - 
and Oetiver - ^ ^ . * ' ' 

We have discussed the experimental design, nature of the sample, defiinltlons *_ 
of variables ; and so forth earlier (Hannan^ Tuma, and Groeneveld, 
1976, 1977). Here We use. similar data that cover a period of 24 months 
rather than 18, We focus on estimating the effects of the three levels 
of'inoome support (or guararitee) levels and the length of the treatment » (3 or S ^ 
"years)* Because the design involved 'stratified, random assignment to treatments^ we* 
include a number oi pre\xperimental variabres in the models: dummy 
variables representing 6 nopnal income categoriej, a site dummy (l=Denver) , ^ 
a dummy variable^ for previous AFDC eKperienci, number of children, a duimny 
variable for having any children uader^ six, the woman-s age, her wage rate^ 
and her years of schooling. ' " 

Though other causal variables, as well as the experimental treatments^ ' 
may have different effects on ^Ates of events in differ€nt time periods, we 



are primarily interested in the time-dependence of creatment effects. ' 

There are a variety of reasons for' expecting the effects of experimental ' 

^ ... ^ ^ ^ A 

treatments to vary jver the experimental period. Reasons for treatment 

effects to be smaller initially than later on include the possibility that 

subjects do nop fully understand the treatments Initially and that.ihey may^ 

need Co search for an opportunity to change state Ce*g., , to find a marital 

partner). Reasons for treatments effects to be larger initially tharf-«ter 
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on Include the possibility that the Creatments become less splient . a^s the 
end of the eKperiment approaches and the possibility that the treatments 
.dmprovfe opportunities for a changfr for some fraction on the verge pi, changing 
marital status- before the experiment. ludirect treatment effecti (e, g.> effects 
of treatmenti on- work ^behavior J which in turn affect marital status) also * 
can cause trEatment effects to vary over timej but it cannot be predicted 
a priori whether ^indirect effects enhance or dampen the initial' response* Our 
inability to predict the shape of the pattern of time variation is, of cours.e5 
the reason for subdividing the obseryation period and letting rates. , ^ 
vary from one period to another. 

In tlio time^dependent models we began by ^aing four time periods with 
end po^tnts of 0.5s 1.0, 1.5, and 2,0 years after the start of the eKperiment* 
These represented a compromise between two conflicting goals. First, we 
wished to have the number of observed events per period be large so that 
the standard errors of parametfrs would be comparatively small, " Second^ ' ' 
we wished to have a large number "of time periods so th4t we could detect 
the shape of the pattern of time-variation (which according to our reasoning 
might decrease over time, rise over timej or riste and then decrease, etc.) 

To improve efficiency of the estimation of effects still further, 
we estimated an equation that allowed treatment ef f ects , but not the effects 
of other, causal variables, to vary over time. This model can be represented 
as follows \ 

jkp ^jk- "jkp-TD ' . p-1 , p ' 

where X is the vector of other causal variables, E is a vector of experi^ 
.mental treatmenc variables, and p refers to oni of p ■periods. 
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i , and each 8., GOn£ain a constant term; however^ for each transition 
^jk -jkp ; 

' - w ■ - ' ■ ' ■ 

j to k| onlyP of these P + 1 constants can be identified* Therefore, *to 

9 



1-dm ntlf i gap±on-^we--mrbi-traTl-ly-^o 
-The P conetants permit^he rate of 'the control group to vary from one time 



period .to another, even though the affects of other causal variables 
* ' 10 

mmy not, ^ . ^ . = 

1.. ■ - • 

The method of maximum likelihood can be uocd to agtimate para^GtGro in 
- . I. ' I' 

equation (43)* The likelihood equation resembles (39) except that there is 

a different expression for the^survivor function and each transition rate has 



an additional subscript p to indicate the time period to wjiich it applies* 
(For <fWails,''s^^ the Appendix.) In other respects procedures for estimafting 
and testing the model and effects of individual variables are the same as for 
a mo^el with transition rates that are constant over time* ^ 

Results of the estimation could , be reported many different ways* 
Because the effect of a variable on the rate itself is of more interest \ 
than the effect on the log of the rate^ 'we report the antilog of estimates of 8s 
The antilog indicates the multiplier of the rate for a unit increase in 
a variable. For dummy variables, which we use to represent experimental 
treatments, the antilog of the coefficient of the variable is the ratio 
of the rate for those whose value on the dummy variable is unity to the 

rate f^ those, in the omitted category, For example, if for some ^ ' * 

6 " ■ 

experimental treatment, e" - 2, where 9 is the coefficient of the dummy 

variable representing,, the experimental treatment, then the rate for those 

on the treatment is twice the rati for those in the control group. The ' . . = 

! ■ ■ 

4% ■ ' 
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percentage change lii the rate for an axperlmental treatment relative to 

■ ' ^ 8 ' 

the rate for those in the control group is Just lOOCe -1>* Thus, 

ft ^ I . .. 

- 2^ ^^ he-pereent^ge-ehange ^ — 

' ^ ' ^ * 1 ■ ' -.^ ' . . 

to the rate for the control grqup ±m 100%? 

Results for T^e-lndependent_Models ■ , . 

Table 1 gives the results for ^'Cime^independent models of rates of marital 

■J- ■ ^ ' * ■ 

status change and attrition. Wa have discussed similar results and their 

interpretation at length elsewhere (Hannari, Tuma and Groenevelds 1977; 

11 ^ 
Tuma, Groeneveld and Hannan^ 1976), Here we comment briefly,, concentratirig on 

the effects of the support levels, . ' 

(. V Table. 1 about here 
All four models significantly (,001 level) improve upon a constant rate 
.model. The set of experimentar treatmanti significantly improves upon a 
model that includes only the other causal variables in the case of the marital 
dissolution rate (*001 level), the attrition rate of married women (, 10 level) 
and the attrition rate of unmarried women (.05 level) s but not in the case of 
the remarriage rate, ' 

Women in each support gro,up have higher marital dissolution rates than 
comparable (in terms of values of other causal variables) women in the- control 
group* Effects of the support levels are large (ranging from a 57 to 129 percent 
increase in the rate relative to comparable controls) and statistically 

significant at the .01. level, except for the highest support* 

■ * ' ■ . 

As expected, attrition rates of those with a financial treatment are 
lower than those of comparable controli^i ex,cept in the case of married women 

" " ----- . J. . - . j^^ ^ _ ^ 

on the low support* However, there is no clear pattern to the effects of the 
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. ^ ■ \ ' ■ 27 " . . 
different support levels on attrition, and none of the Individual cpefficients 

/ Is significant. ■■ , ' . * , . 

Reiults for the Time-dependent Models ' ^ 

We have estimated time-dependent (four-period) models of rates\of 
5^ dissolution^ remarriage and attrition. See Tuma, Hannan and Groeneveld (1977) 
for a more detailed dJLscussion* The results of these analyses indicated that 
eKperimental effects oh attrition rates did not significantly vary over' time < 
or in- any patterned way, Experiment^l= effects on the remarriage rate were 
significant at the *10 level, but there was no particular pattern to these 
effects, suggesting that they resulted fro^i chance alone* On the other hand, 
experimental effects on the disaolution rate had'a striking pattern. So : 
we focus on these. 

Table 2 shows results on the experimental effects on the dissolution 
rate when the total observation period (24 months) is treated as a single 
period, four six-mont-h periods, and two periods (the first six months and * 
the remaining eighteen months) = In all three models the set of experimental 
treatments significantly improves upon a model that contains only the other 
causal variables . 

T^le 2 about here ' , ■ 

-J , ■ . ■ , '' 

Support level effects over the four periods are plotted in Figure 1. The plot 

■ shows that all ^ three support levels produce an exceptionally large increase 

the marital dissolution rate during the first half-year of the experiment,^ 

Except for the $5600 support level in the second^ half-year ^ effects of the 

support levels in periods two through four are positive' and show no clear 

;= pattern of variation over time. This suggests that (1) the support levels 
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have traniitQty effects on dissdlutipn rates of white women that subside 

after six monthi^ b^ut (2) they also have nontransitory positive effects on 

the dissolution -rate, furthermore, the effects of each support level in che^ 

first six-month period^relative to its effect in each subsequent perio<^'?s^ 
• " • ^ . • ^ ■ ' 1^ " ' — ' 

about the same for all three suppSrt levels, ATthough there is a clear ' 

' i . 

^ ^pattern of time-varying effec^ts, the. four-period model does not significantly 
improve upon the one--period model. 



' / Figure 1 ^bout here . ' 

^ With these findings in mind, we estimated a two-period model in which - . 
the first half --year is distinguished from the rest ,of the experimental period, 
Rel^t^ve races across support levels- were constrained to t^e equal in the 
two periods (see .Table 2), but their rate relative to controls was allowed 
^ to vary from one period to the other. Treatment effects in this model 
significantly (.05 level) vary over time, according tq both the likelihood 
ratio test and the F test on a dummy variable for an effect of financial = 
treatments in the first half-year. 

The estimates for the two-period model (Table 2) indicate that the 
. effect of each support level is 2.32 times as aarge in the first half- 
year m it is thereafter. Based on our prior^/a^uments , , this suggest 
thac the treatments immediately changed the oji;^brtunities of some respondents 
of the verge of dissolving their marriage. Byt the treatments arso seem 
to have changed the long-range opportunity scructure. The effects of tfta 
support level's during the 0.5 to '2.0 year pSriod are potitive and sig- 



1 .4, 
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nif leant (except for the $5600 support), |hough 14 to 15 percent lower 
^ ^ , 31 ' 



than in the results for tlie tiiiB-lndeptndlnt model, Sind^^the long-range 
impact of 'an InQome? maintenance program should not depend on transitory 
effects on rates during an Initial adjustment period ^ the eff feces of the 
upp©^fe^^ve-l^--in-t^he-Q>^-t 

estimates of the ultimate effects of income maintenance than estimates based 

if 

on the one^period model. 

How Well Does th^ Model Fit? ' . 

Markovj^models have a reputation for fitting data poorly* But our 

' - * - 

extensions of the Markov model should have h^elped to improve the ability 

of the model to fit the data. : So' far. we have used 14i^.elihood ratio tests . ^ 

to assess the relative fit of a series of nestM models. We havi' learned 

that some models do not improve upon others. Here we look at the 

absolute fit of the "model and compare it to a common alternative I 

Three main questions are involved in assessing our dynamic model: 

(1) To what extent do predictions based on our model differ systematically 

from observed values? (2) If there are systematic diff erences , are these 

related to the experimental treatments, and so causing us to make 

erroneous inferences about experimental effects on a particular outcome? 

(3) As compared to approaches that seek to minimize prediction errors 

for a single observable variable, how well does our model explain sample 

variation in outcomes related ^to marital stability? 

In answering these questions we consider three kinds of outcomes* 
the probability of being in a given state (e.g.j single) at ai;iy moment, the 
expected number of marriages and marital dissolutions in any given time 
interval, and the probability of leaving the original marital status in 
some time period. Our model implies that »each of these is a function of 
marriage, disssolution, and attrition rates, as indicated in Section I. 



'Thus we ean use a woman's values of treatmenta and other ^cauial, variables • 
and the e^timatad effepts of these variables to predict her rates of marrl- 
«age, marital dissolution^ and at'trition, and then her values of "the^ differ- 
ent outeomes listed above. We chose to consider these outcomes a^t two 
arbitrary^ times I one and two years after the start of the experiment. 
For each woman used in our analysis of trahsition tates^ we retrieved her 
' Ob sarvafl^ values on the outcomes liited in Table 3, : These variables ware not 
directly used in estimati|ig transition ^ates; of course^ they are not independent 
of those data either^ which is the reason for having a single model. We also 
predicted' each woman's values on these outcomes using the estimated effects 
from the one--period models of remarriage and attrition and the two-period model 
of marital dissolution (see Tables 1 arid 2). ^ ^ . . 

Table 3 about here 

To detect systematic differences between predictions and observed values 
we report the mean residual for each outcomes i.e*/ the mean difference between 
observed and predicted variables. We also report the observed mean of each 
outcome because the relative size of a systeiriktic difference tu^.ot some interest 
too* With predictions from a linear regression models the mean residual would 
be zero. This need not be the case with predictions from our model. The re- 
sults in Table 3 show that the mean residuals for our predictions are usually 

■is. ' I " , 

small bpth in absoluta terms (the largest is ,02) and in relative terms. There • 
is little indication of any overalj pattern to these differences, except for 
the last four; dummy variables, which have consistently positive mean residuals. 
It is well-known that^no change in status has tended to be underpredicted in . 
sociological ^applications^ of Markov models (e.g., see Blumen et al,, 1955), 
Our introduction of population heterogeneity has made this a, comparatively 
small problem, but it has not erased it entirely. 



Small mean^residuals could Mdi lyitematic differeneea aisgeiated with 
different treatments I whlcfi is clearly undesirable if con^Toi-f inancial 
differences on an outcome are of interest. To answer the second question, we 



peffor:med ons-^way ANQVA on the rasiduals for each outcome* Treatment 
dlfflirences in the residuals never even approached statistical .significance , 
(the smallest prob-value was^ greater than \50O 

In addressing the third question, we focus on a single, coiranon inexpensive 
alternative----linear regression analysis* We regressed each observed ontcome ^ 
on the prediction from our model; we also regressed each of the fourteen outcomes 
on initial marital status, treatments, and other causal variables used in our 

analysis of the^transition rates. For both our modal and the regression model 

2 ' 
we report R*"* the square of the correlation between the observed and' predicted 

yariables. Since we expected a poorer fit from our model than from one designed 

to minimize errors, we were surprised to find that for ten of the fourteen 

outcomes our model explains more of the sample variation than does lineiir 

,- ' -'^ 
regression analysis. Moreover, the advantage of our model in these ten cases 

tinds to be larger than the advantage of linear regression in the other four. 

^ Though we have cd^sidered tha predictions of our model for several outcomes ' 
at, two arbit-rary times,, we have not yet seen how well it predicts the time 
path of these outcomes* Computajtional expense has forced us to 'examine the time 
, path of only one outcome. We selected /ithe proportion who are unmarried at time 
t, conditionaL- on having not attrited by time t. We chose this because it is 
similar to the most important policy outcome, because it should reveal whether 
' "experimental effects are confined to an initial, brief adjustment period," and 
because it depends about equally on our estimates of marriage, dissolution, and 
attrition rates* This choice provides a severe test of our method because 
it yses estimates of all four rates of change (r r^^, r , fo^)* , h 

i^ iJ 21 ^> 



Figure 2 gives observed and predicted curves for this outcome\by 

support level. Points on the observed curve are given by N^j (t)/(l-N2j (t)) 

where N (t) and (t) are the number in treatment j at time t who are' 
IJ 3 J * 

niBawied-ftnd-mttr4ted-5-~re#^ 
are calculated as , . ' = 

J 1-1 - - , 

where N ^ N (0) + . (0) (tho initial number in treatment j).^ and p ..(t)^ 
V j Ij 23 ^ ^ iji 

and , . (t) are calculated for each woman i enrolled on treatment j using 
- 3ji - 

equations (28) and (30)* respectively. Predictions' ar,e based on the on^- 

period models of remarriage and attrition and on the, two=period model of 

dissolution; they assume eaqh woman has 'her assigned treatment. 

^ ■ 
Figura 2 about here , 

We 'begin by considering the observed Csquiggly) curves. ' Fitstj note 
that the proportion of unmarried women at the start of the eKperiment 
differs greatly from one treatment to another. These initial differences 
result from the use of a stratified random design in which marital status 
was a •^stratif Ration variable. Unmarried ,wotnen were more likely to be 
assigned to treatments with a lower support level. Because of this, a 
comparision of ''posf-^test" levels is clearly inapproptiate. 

Next, uotice that the observed curve for the control group is relatively 
flats suggesting that there are no important "natural" time trends (due 
to agings secular change 5 etc.)* On the other hand, the observed aurves 
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for the' three support levels show ndticeable increases in yie" conditibnal 

proportion of unmarried ^omen aftet? two years (+*039 for the low support i 

. z:..^ . ■ r__j^ _ _ - . ■ ■ . : • ..... ^. i . ' '_ . ' 

+,044 for the medium support^ and +»031 foe the'hlgh support versus 

+.009 for the cdtitrol group). Furthermore, the p^roportlon of unmarried 

women among those on the financial treatments rises fairly steadily 

= throughout the two--year period.* This upward trend is, quite apparent for 

the dow and medium supports. It is less certain for the high support^ which has 

the fewest subjects and the most extreme fluctuations about any overall 

time trend. There is little evidence that the proportion of unmarried women 

among those oti financial treatments has reached a plateau within the first two 

years of the experiment , as we would expect if an equilibrium was reached during 

- this period. This suggests chat dynamic analysis is really needed to assess. 

tKperimental effects accurately. We elaborate on this in , Section III. 

. Now let us consider the fit between the (smooth) curves- predicted by 
our model and- the actual (squiggly) curves. We rely on visual inspection . - 
to compare the *two sets of curves. On the wholej^^the fit is quite good 
■ except for the high support groups for whigfi^ the actual curve is noticeably ^ ^ 

below . the predicted curve. Because only 240 women are in ^the high support 
group, a change in status of a very few women makes a subatantiai difference 
in the observed curve* . Hence phe deviations for this grQ^^Mr^.- less 
worrisome than they would be for a large group like the controls . (N-847) , 

^ . - - . ■ 

Our scrutiny of the implications of our model for various dutcomes 

at arbitrary times has revealed no major disadvantagjs and even some 

small advantages. The model primary advantage is^ of course, it^ 

ability to predict the ^ime path of a variety of interdependent outcbmes 

^ reasonably well, 

ERIC . ■ 3e . 



Ill COMPARISON WITH OTHER APPROACHES - . , 

' We return finally to issue raised at the outset i the advantages of 
event-history analysiai relative to other ^approachesS We continue to assume * . 
throughout this discussion that events are gene^rated hy a Markov process . ^ 
whose transition rates are log-^linear fupctions of f^ogenous variables (36) , 

dbviously any comparison of alternative approaches depends on the 

-- . . . - . " ■ . ^ 

assumptions about the pv'ocess generating events* Following Colfemaji's 
(1964 J 1968) work, nonetheless, we wish to challenge thdS itill 
widespread view that, substantive assumptions ought to be dictated by the - 
form in which, the data are collected. Th|it is, we do not believe that- 
sociologists ought to change their assumptions about the lender lying process 

(their model, in our terminology) when they shift from analysing panel data j 
"sayj to analyzing cross-sec^tions or ev#nt histories. In our view^ a major 
advantage of formulatdng probleias In terms of dynamic stochasxic models is 
that we can use different data structures to estimates parameters of the 
^same model. This provides a way of unifying a variety of data analytic 
procedures. : . . . : - ^ . ^ 

We begin with' ^n extended discussion of ^cross-sectional ailalysis because 
this has been—and will undoubtedly continue to be—the mainstay of sociological 
rese&rch. .We then contrast event-'history analysis with two other strategies 
that use .temporal data: event-'Count analysis and panel atialysif ,- ^Our. 
discussion of each is brief. To the best of ^ our knowledge, event-=count 
analysis has* not yet 'been developed (let alone applied) except for the most " 
elementary kind of Markov model (a Poison model> . So^ our cc^^ntp ori,lt *fre 
intended to encourage the developmentf of this approaoh. On the ofWir;handj 



Singer and. Spilerm^, (19745 1976) have treated in detail the difficulties of 
partel analysii witK categorical dependent variables.^ We mainly review theie 
^to emphasize that they -are either absent or much less seripus with eyant- ^" 
h^^ry analysis, - . 

Cross^ieGttonal .^alysls - • ' ^ * , 

'Grbai^sectional data give the sGate that each member of a sample . 
^ccupiei at a particula^ time t* Earlier we referred to the unconditional 
probability of being in a state J at time t as the state probabilitys p,(t)* i 
Given n possible states, there are only (n-^1) unique state probabilities ; 
since the *n probabilities must sum to^unity* 'Butj in general, there are 
n(n^l) unique transition ratps. Bacause ft(n»l) > (n-1) for n > 1, it is 

ipmediately obvious that cross^sectional analysis does not allow altt parameters 

. ' " ■- 2 ' ' ' 

of a model to be identified uniess we can specify (n^l) of the transition' 

rates, elther^on theorbcical ^rgy^ds or^rp^ a priori knowledge. With event- 

history analysis 'we ^ can estimate all parameters. So eVent^history analysis 

is clearly preferable to cross-'septional analysis if we wish to understand 

the process fully or to predict other outcomes. 



Under certain circumst^ncesV however, cross=sectional analysis can si 
useful "information about the process gehirating events^ It is worthwhile 
identifying these conditions. We begin by considering the situation in which 

the process has been operating a comparativaiy long tiiAvso that the distribution 

. ■ ■ ^ .... " ';■ V : . . ^ 

of the population across states is in equilibrium. For concreteness we' again. 
^sCart, with : the twp-^state ^case,' ^ - v 

Since the rates^ r^^ and r^^^^s caririot be negative, the two-'State model , - 
implies that the^ probabilities of being in each state eventually .reach stable, 
so--called stQady^state, values: ' * * . - ^ 



. .^12-21 %2\ =21= : ^ 



2 ^12+^21 • 42^, ^21^ ' ~ ^ 



Of course, indivlduais continue to ^anga from ohe state to the other. That 
iSs the model implies an equilibrium ^probability distribution on the ,a^gregat% 

. . • > / " ■ " • . •^v^.Jv.^-, ■ 

level; it neither assumes nor impliii equili^^^um^ on the individual, revei*^^ ' / 
\Note that.^ . . w . ^ ^ . \ 
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where " 5-2])— 12' exaj^Jle, this^model implies that in the steady-stata, 

the log of the odds of being . dfnmarrie^i" ^rathpr than married)^ is lineat ii\ X. ' 

Equation (48) is the usual form of a binary logit model (Berljsonj 1944; 

^ \ . \^ 

Theil, 1969j 1970), and if all members of X are dummy variiMeSj then it 
is just a special case ^^f Goodman's (1972) log-linear model of the odds 
ratio. . Thus> when a population is'"ift-. equilibrium, logit (or log^linear) ..p^ 
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analysis of ■ cross-sectional data tells us the difference in the. effects of 



variabl'e's *oh^the. ttwo-^ Note that "no effect" o£. a variabia 

ih' thi crbss^^'sactlbh^l ■ ibglt analysis can be due to iti equal effec'tp dtf the 
rates .^i^^hajt^s^t^^ 



ifrevelant to the pro;ceiSs only t^at it has no net ef feet oti the steady-^ 
states disitibutlx^n, ^ ' / . ^ ? ^ / ; . . 

^Ss there a aiftilAi eonnaction between the general^ n--state ^torkov 
model and ^ttultWpmial Ibg^t analysis? Unfortunately, the answer is n,o* This 
can' b'e proved by a, ^single contrary tassi Consider In (^) /p« C^) ' for the' 
/thfee--state mpd©l w^^ attritloni^'^hich w.e have used in out empirical 
^^talysiS, ' though both p (^) and'^piC^) are zero (because "'Eventually evetyone- 
"attrites^") 5 they do have a f inite ;tfatio * ^ We find that^i'^N.. ' ! ■ ; " 



The ^jl^ 



ssion 0^ the^ right-hand side ' if quite cpmplicated afid substitutlpji ! 
of eqt|ations like .(36) .daes iiot produce arty thing' resembling -a logic modei/.in 
gener-fil.^' . (If . r^j^t^^ i-e.^ attrition rates for married and unmarried vW'b^ien 
are iden|ttcalj then (49) does simplify. to (48))^ 



However s an impgj^^ant class of Ilarkov models^^ feniiral birth and death 
mod ^1 s "^df hav e a t ^^^^K" s ta t;e_ d is &r ib u t ion that has ' t he^f o rm . of ■ ^ . 



multinomial i^^^^^o'He|<r In these models spates can b^ ordered; that 
there, can oi^&r^'only transit^ions between neighboring states .(e.g^^ in a 
/three^state model transitions from 3 and from 3 to 1 are impossible). 

Then, for examplpV in a three-state model with ordered states, equations. (47) 
and (48) continue to apply ^ and in ad^tioni = .,. 
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(50) 
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^^;;.:l«f-(Ey-^:-^23i ' , : . ■ (51) , 




■ f t52) 

The similarity, in form of (52) to (47) and (50) maans that there are no 
clutea In..^ croas-sectional data to tiXl ^us how to order tha states. This 
must be dona dn theoretical grounds— or else we must have eyent-^his tory - 
'datas which do permit us to obserYi^ what kinds of transitions cam occur. 
But, if we do .know the order of^taces in a general birth and death model, cross^ 
seGtional multi^#omial logit analysis does let us make conclusions^^ ^bout the net 
effect of a variable on transttion rates between states. 

. ^ So far we ^ave considered .the situation in which the $ystem is in * 

; n ' ■ , . , 

equilibrium* We have indicated that in the ^3teady-state, the- log of the v_ 

fi' , ^ 
odds of beiftg in one state rather than another has'^-a very simple form if 

states can be ordere^. On thfe dthar' handj if the steady-state has not 

been reached, then the log of the odds of being in "pne state rather than 

another is a very complicated function of time, the transition rates ^ and 

the initial conditions^^even for the very simple two-state case. This 

p ' ' ' 

can be seen by forming the ratio of the right-hand sides of equations (26) » 

and (27). 



■'Suppose we perform cross-sectidnal* analysis at two successive time 
points and determine that variables have very different effects, We 
cannot be certain what to cpnclude. Has the; underlying process changed—that 
iSj'has thd relationship between variables and transition rates altered? 9 
Or has the system just moved closer to its ultimate steady-state value? 
Without some form of temporal analysis we cannot answer these questions. 

Social scientists have been so wedded to crois-^sectional analysis that 
they seldom seem to reflect on the likelihood that equilibrium exists or 
on the length of time required for a system to reach -a 'new /e'quilibriu^ 
following soma intervention or structural upheaval. "'/We suspect that in 
most ca'ses inertia greatly slows the speed with which social^ systems reach 
■equilibria. We note that the equilibrium assumption implicit in sociological 

theories prominant a few decades ago (e.g.^ f unctionalism) began 

to ,be attacked more than a decade ago. But these criticisms have barely 

begun to penetrate sociological ■methodology* 

We will make our discus?ion of this issue more concrete by referring 
to our empirl^ckl'^illustratlon . Many analysts faced with dWta. 
like ours\.,:pigjit ^conduct some sor.t of logit analysis that ass^amfes that 
equilibriunf' is'" reached within the experimental period. But can we expect 
the steady^atate probability discribution of marital stitus to be approached 
during the 3 or 5 years of the experimental period? As social experiments 
go, SI>ffi/DI>ffi is long 5 so it might seem that this would happen. However, 
according to the models w^have discussed, how Idngit^ takes to approach a 
new steady-estate (say, to within .01 of itiv ultimate value) depends on both 

preexperimental and experimental rates of marital status change. Our results 

« " " ii - • . 

imply 'that SI>ffi/DI^ is much too short for the steady-stace to bj^- approached 
during the experimental period. 
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■' .Marriage and dissolution rates' are apptosilmataly equal to annual 
probabilities of fbrmlng and dissolving a'^marriage, respectively, Among 
•SIIffi/Dii^'a participants, who Have incomes below the U.S. median, both rates 
are sdrafewhat higher than in the overall U,S, population. In the environment 
facing the coritrol. group, - ,10 and r^^ * ,05 are fairly typxcalJ^ 

If an aggregate-level equilibrium exists at the beginning of the experiment, 
. the initial probability of being unmarried is .05/(,lp + ,05f ^ ,333, 'This 
. is close to the observed proportion of unmarried women among the control 
group at the beginMng of SM/DIM (see Figure 2) , ' 

Our analyses reported in Section II indicate that for whites SIME/DIME 
has a negligible effect on marriage rates, but ^roughly doubles the dissolution 
rate of chose on most treatments. If the dissolution rate under incoma, 
maintenance is twice that of the controls, i,e., r^^ ^ *10, then according to# 
the two-stkce model the equilibrium probability of being unmarried under 
income maintenance is ,10/(,10 + .10) ^ .50, Thus, the model predicts that 
under such conditions j the proportion of unmarried women in the population 
would eventually increase about 50% above its pre^income maintenance 
value (from .333 to .500). However, the proportion of unmarried women 
would only increase by about .04, or about 13%, in the first two years 
(see Figure 2^), and by about ,09, or about 26%, in the first five years.: It 
would take nearly 19 years to be within .01 of the steady^state proportion. If 
data on marital status of participants at any point during the 3 to 5 years of 
the experiment are analyzed cross--sec tionally , the ultimate effect of income 
maintenance will be greatly underestimated. As we mentioned earlier, the 
observed curves for support lavels in Figure 2 do not suggest that a plateau 
or equilibrium has been reached within the first two years. . ^ 
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The only way to decide whether a system is in equilibrium is to 
collect data over time and to analyze it dynamically—that is, in a way 
that lata us study the timo path of chan.^e in the phenomenon. 

Event-count Analysis 

There have been a number of sociological studies that analyze the number of 

times a particular event occurs in some time period. For example, Spileman (1970^ 

analyzed the number of riots per city during the mid'-sixties • We refer to 

this type of analysis as event^count analysis* We suspect that a number of 

surveys have asked such questions as "How many times have you been married?'* 

and "How many t^lmes have you been divorced?" However, we are not aware of 

any event-count analyses in sociology where counts of more than one type of event 

are analyzed within a single model. 

*• ' 

Given the assumptions of a Markov model we can derive expressions for^ 
the expected number of different types of events in ^qme time interval* For 
example, for the two- and fehree-state models discussed in'^this paper, we 
derived the expressions given in equations (31). through (35) . These 
equations, combined with observed data on the counts of different types of 
events, permit transition rates Co be estimated by a nonlinear regression 
approach. That is, we can estimate rates^ — or the causal effects of variables 
on rates==-by minimizing the sum over all units of the squared deviation between 
the observed count of events for each unit and that predicted by the model. 
This approachs which we have not yet used', has one inherent limitations we 
know of no theorem (comparable to the Gauss-Markov theorem in linear 
regression analysis) that estimators obtained in this way will have optimal 

4 

Statistical properties^ — even in an infinite sample. 
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^Maximum likelihood estimators typically have optimal asymptotic proparties.. 
But to perform maximum likelihood estimation we must know the probability mass 

mts. To the best of our knowledge, the expression 
for, this furiction has not yet oeen derived for a general n^state Markov model. 
In fact, it is not even clear that such an expression can be written in closed 
form for a general n-state model. The probability mass function for the number 
of events can^ of course^ be written explicitly for certain special cases (e.g., a 
Poisson mouel) , but it is mathematically intractable even for the two- and three- 
state models used in the illustrations in this paper. 

Given these difficulties, it seems obvious that event^history analysis is 
preferable to event-count analysis. Nevertheless, event-count analysis deserves 
further study. Under some circumstances, event counts either already exist or 
are feasible to collect, while event histories or panel data cannot be obtained* 

Panel Analysis 

■ Panel data, which record the states occupied by members of a sample at 
a series of discrete points in time, are the temporal data tria^t commonly 
available to sociologists. Singer and Spilerman (1974; 1976) have identified 
the following problems regarding estimation of transition rates in a general 

First, sometimes panel data on categorical^ variabi cannot be embedded 
in~that is, described by—a Markov process. Moreover, sampling variability 
and measurement error' can cause panel data to be unembeddable even though they 
are truly generated by a Markov process. Sacond, even though the panel 
data may be describable by a Markov process, there may not' be a unique matrix 

/ 
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of transition rates that describe the data* Furthermprej the different 
matrices obtained in the nonunique cases may suggest substantially different 
qualitative conclusions. (See, for example , Singer and Spilermans 1976s p» 31) 
Neither embeddability nor uniqueness is a problem in event-^history analysis 
because maximum likelihood estimators based on such data give unique 

o 

estimates of rates-=or of causal effects on rates. 

Third, Smger and Spilerroan (1976, pp. 44-48) note that small changes 
in an obsery^ed matrix of transition probabilities ' (due to sampling 
variability or measurement error) can sometimes ^lead to very marked 
changes in estimates of transiti'on rates* This is clearly undesirable* 
On the other hand, in our experience in analysing event histories , giv^n a 
moderately large sample, fairly substantial errors- in records on the 
occurence or timing of events do not qualitatively alter estimated patterns 
of causal effects of variables on transition ratesp Insensitivity to 
sampling and me^^rement erro.r is, ^we beiieveV 4fi ^important advantage of 
event'-his tory analysis. Because such eryors ^re' unavoidable ^ the sensitivity 
issue clearly deserves further study—in both panel and event'-history * 
analysis, , ^ * 

Fourthj estimation of transition rates from panel data is also sensitive 

/ 

to the length of the time interval between waves of the panel. When the 
time interval is large, each row of the matrix of transition probabilities 
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approaches the steady-state probability distribution. (See-, for eKample, 
equations (8) through (12)). In this situation there are only (n«l) unique 

= transition probabilities, rather than n(n-»l) * This means that the data 
contain no more information than cross-'sectional data. On the other hand, 

•if, the time interval between waves of the panel is very short, almost all' 
members of the sample will be in their original state. This is not very 
inf orinative , either* With event-history analysis the length of the observation 
period cannpt be too long. It can be too short—if no event has oocurred. 
However, Tuma and Hannan (1978) have shown that with samples that are moderate 
in size, rates can be estimated well if as few as ten percent of the sample 
have had an event. We have not seen comparable results based on panel analysis , 
but we suspect that it does not perform as well. ^ 

Fifth, we know of no way of estimating parameters in; a general n--stace 
Markov model from panel data when transition rates are functions of exogenous 
variables, as Is almost always the case in problems of interest to sociolo= 
gists. Singer and SpilarTnan (1974) have reported some work on estimating 
paramecers from panel data generated by a mixture of Markov processes, when the 
mixture is described by some specified probability distribution. This work 
is helpful, but it still doesS^t permit causal inferences to be made. As 
we have shown, causal relationships are easily studied with event^history 
analyiis , 

A sixth problem with panel analysis concerns the ability to study 
and detect time-dependence in the process generating events. This appears 
^ to present a very difficult problem for panel analysis. As we have indicated, 
various kinds of time^dependence can readily be^invejtigated through avent- 
history analysis. 
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Conclusions ' ■ . 

Our CQncluslons are very simple. Event histories provideV,E'ich 

opportunities fo^r answering fundamental sociological questions. We have 

shown how to anaiyze event histories when data are generated by a well-^^ 

behaved stochastic^' process , ' The procedures we have outlined permit analysis 
} 

of causal effects on the rates at which events occur and_^f .time -dependence 
in such ifat^,' These grocedures are simple to implement, and in our empirical 
application they yield good predictions about a variety of observable variables, 

^ Event=h±story analysis offers substantial advantases over other common 
approachas to 'the study of causal effaces' on changes in' qualitative variables'. 
Since ,4n many situations it is no more difficult to obtain Information 
on cte timing of events than Che count of events, we urge that sociologists 
begin to' c'*116'ct and analyze such data. r-" 

\ ■■ . ' 

'■-^ ' „ ■ 
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■ senting the opinions or policies of the states of Washington or Colorado or 
any agency of the Unitted States Government* 

2 

A simultaneous- linear , equations model might appear to be one alternative* 
However, the same exogenous variables appear in every equation, leading to 
underidentif lea tion of the structiirai parameters. 

^For a more thorough discussion, see CoK^and Miller, 1965; Feller, 1968; 

i \' 

Karlin and Taylor, 1975; or any other standard text on stochastic processes. 
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, ■ . . . \ ^ • 

Mathematically these assumptions can be expressed as: 

0 1 Pjij(ii.t) < 1 ; • - . 

P(t,t) - I . ' , V 

5 , 
Equations (8) through (12) do not apply to a general three-state modml, 

i^e,, one in which r^^ and r^^ greater than zero, i^^^. 

6 

As m&ntioned by an anonymous referee, we could have chosen to model 
r and m., t.,lt.) rather than r., as functions of exogenous variables. 
The choice involves a decision about the substantive nature of the process; 
it is not Just a methodological issue. The 'Approach sXiggested by the 
referee, which has also been advocated by Tuma (1976) and by Singer and 
Spilerman (1974), J.s appropriate when the decision to leave the current 
state and the choice of a destination are separate or seq\ential. In the 
example used In this paper , it seems reasonable to think that the decision 
to leave' a state and the choice of a destination are no t separate; therefore, 
we have not adopted this other approach. For other social processes it may, 

indeed s be preferable to model and m separately. The method of 

. . . y 

estimation and aGcompanying soft-ware described below can readily be used 

' ' ~ • ■ w ' 

to estimate r^ = from data on the dates of entering and leaving state j * 

^ A FORTRAN computer program called RATE has been developed to find the 
maximum likelihood estimates of parameters in Equations (35) , (37) ^ and 
(43), among otherd. Written documentation, test data and test output are 
also available, For information^ write the first author. 



For our definition of marriage and marital dissolution, see'Hannan, Tuma, 
and Groanavfeld, 1977. ' jf 

The same problem arisfs' in regression analysis when dummy j^^tiabltSi/' 
are used to represen^t a categorical variable* One categoryfemusCv^B^^ 
omitted to achieve/identif ication; . the one chosen affects the interpre- 
tation of coefficients of the included dui^y variables but does not ' 
affect predictions of thfe dependent variable^) 

10- - ' ' ' 

The rate^ of the control group may vary oyer time beGause of aging, secular 

trends , .fete . , ' ' ' ' 

11 - - 

Althoifgh these papers do not explicity mention attrition as a third 
staCe,: the* actual estimation procedures were the same ,ms those reported 
here , 

■If there is no attrition, the probability of a change in.marital status' 
before t is 1 - e:<p[-r,, t]. (See equation (17))* ^ a Taylor. eKpansion 
this is 'approximately 1 - [1 - ^j^^l = ^jk"" J^en t " 1, the probability 

of a change is approximately r , ' , 

13 ' r 

These numbers are obtained by rounding off the crude prop^ortion of co^^rols 
who marry (if initially unmarried) or end their marriage (if iti^^ially 
married) in the first year of the experiment. ' ^ ^ 

14 

Some of these problems do not arise in car tain special cases , For example, 
the first, second and fifth problems mentioned below do not occur for the 
two^state model. ' ^ / 



We used this' feature in Hannan,' Tumai and Groeneveld (1976) to eliminate 

' ■. 
the possibility that effects of income maincenance Creittoent on the' rgite 

marital dissolution were due to attrition bias. Our unpublished Monte 

Carlo studies show that random error in the timing of events has surprise 

ingly small effects on. the' 'qiiality of maximuni likelihood estimators of 

ra^fes qbtained from event history analysis-.;^" ' . - 

- ' ; ' ■ ■ 5 ■ 

■ i ' ■ » -3' . » , * » ' ■ 
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, All equaciuns cQntaln^i^ther cauial vaflabka gfii|hlnTabli L Coiffklenti are explSj'^p IndldiN the Qultlpliefi'of *£hi rati., A' cwf f idint 
. ofli fe^iiiiH tiio iiliibli? lias no ufruCl 



* OaO ^ 5:5.; ^*{};05 ^ p > 0.0^. ■j;**^!)! 
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i' Table 3, Obocrvcd and Pvodictcd Vnl-joa of ArblCfnfy Outconcn 
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Misn 


Mean 
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1 Q 1 7 *t 
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i- 
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0 * Hirried ac t^l 

■ ittrlteiat t^ • 
Attrltid at t^2 ^ 
^ , Marltai PissolutionS|t^l 
No, Marital Dlssolutiorii,t^2 
Ntf. Marrfeges, t^l 
I p, Mairlages, t«2 
♦Coflilnuously Single to t^l 
ContlnuQusl)^ Single' to t-2 
Continuously Harried Co t'^l 

; f I. 

iContinuously Harried to ts2 



t Except for the variables on mim nf marital uvtMiiN, the oi.oervcd varial)lcfl arc duniniy (0-1) variable 
I There ^ere ?05 Initially single white women' and 1212 initially married white women, giving a total of 



1917 white wonien, 
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FIGUHE i RELATIONSHIP OF INCOME MAiNTENANCE SUPPORT LEVELS TO MARITAL OlSSpLUTION RATES OF^WHITE WOMEN OVER Tlft^E 



^ *0.10 ^ p ^ 0.flii; ' 0.05 ^ p > 0,01; 0<Q1 |)|. 

ERIC . ■ 



60 ■ 



0.60 



0,45 








2 








S 1 0.40 






V. 








_• ' I*. 






dl3S 





0.26 



-CONTROLS 




0.6 



1.0 

TIME I IN YEaI^S 



1.6 



2.0 



FIGURE 2 OBSERVED AND PREDICTED PROPORTION OF WHITE WOMEN WHO ARE SINGLE 

AT TIME CONDITIONAL ON NOT HAVING ATTRITED BY TIME t BY SUPPORT LEVEL 
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ttt^^wi^ertfviikfeiiho^dd equation like (39) for the model with- transition 
ratftg ^i^hat% Yayy p to another (43) , we ne ed the survlyor - 

funct^i^ii; .fqv'' ,statX j^^^ u,X,Z), where u ^ t , (the time of the 

(ra-l):V;;j^ia;ir6itj'^4n^^^ iater times u £ t. ' The survivor funccion 

la oby^ined-by analogue of (21)*, 



% '/ ■ • . 

'j'-'^::' ; ' 'v- :vf;;(t|u,x,z ) ^ G.(t|u,x,z ) 



(AlO 



wher^^.r,;.^'*^:;^ r^Y^^^^^^ )/dt, this implies 



7, 



r.^ dt. 



^Jn^t^graiting both sides from ;u tp t ^ we . o^^ln:* 



':'V^r; 



dt 



cm 



(A3) 



(A4) 



fo.'ielimi.nate the Int 

4eline- u ' and t' a 9 
. P P ^ 



egral on the rigVt^hand' sid,e:-af^-:th^^ ^Kpr^as^/pfi/fl^Dve , 
follows : ■ \'\^y }] 

. ... ■ . "f '" >i 
■ if T < u Of t < T ',•-„;( staf*' j is entered after 

• period p .ends ftr '% -ocl^iii!:!:. b.ef ore fietia'd'^p" Jjegins) 
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u if t . <_ u < t (state, 4 1b entered in period p) 



"^-^X u < T t (state j is entered before period 



p begins ' and t occurs after period p begins) 



t' ^ 0 if T < u^- or t < T , (state j is entered aftjer 

"p. " P ; . - p-1 ^ : " " - > . X ' 

period p Jfends or It occurs before period p begins) 



T if u < T <v>t f(|tate i -is entered before period 



p ends and t occurs after perioa p ^nds) 



t if T T < t < t.'' (t occurs in period p)/i 



where T is the end of period p,^ ^s' defined in-^ the text. Ufing these definition^. 



(AS) 



So, anaidgous 'to (39) , we can write the likelihood equation for the model 
given by (43) as: ^ . . ^ ; " 



^ n N ^ 

^ n n n 

j^l i^l m^l 



G^Ct t ..X,£ )■ 



n 

* - n 



r,, (Xj ) 



(A6) 
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